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Abstract. We use the irradiated secondary star in the 
remnant of Nova Cyg 1975 as a probe of the cooling white 
dwarf. At superior conjunction the flux in the B band is 
dominated by the irradiated face of the secondary star. 
The heated face produces an orbital modulation which is 
dependent on the strength of the irradiating source. We 
demonstrate that the cooling rate of the white dwarf is 
consistent with the theoretical model of Prialnik (1986). 
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1. Introduction 

In many close binary systems the hotter star will heat one 
face of the cooler star. As the orbital motion of the binary 
brings this face into view the observed flux from the bi- 
nary will increase, only to fall again as it rotates out of 
view. After a nova explosion, the hot white dwarf is an 
obvious candidate for heating its cool companion. Proba- 
bly the best evidence is in V1500 Cyg (Nova Cyg 1975), 
where Schmidt et al. (1995) show that the secondary star 
dominates the photometric modulation. They fitted HST 
spectra with a red star whose unperturbed temperature 
was ~ 3000^ but whose face towards the white dwarf 
was ~ 8000.fi:. 

There are three other novae for which there is photo- 
metric evidence for heating. The orbital modulation in DN 
Gem (Nova Gem 1912) was found by Retter, Leibowitz & 
Naylor (1999) to be well described by a heating model 
and DcYoung & Schmidt (1994) suggested heating could 
explain the lightcurve of V1974 Cyg (Nova Cyg 1992). Fi- 
nally, Somers, Mukai & Naylor (1996) found the infrared 
lightcurve of WY Sge (Nova Sge 1783) required a heated 
face to be modelled successfully. However, in this case the 
level of irradiation was so low that it could have been sup- 
plied by the accretion luminosity, as occurs in the dwarf 
nova IP Peg during outburst (Webb et al 1999). 
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This led us to ask if there were further evidence avail- 
able which would help us identify the source of irradiation 
in old novae as the white dwarf. It is particularly impor- 
tant to do so, as the irradiated surface of the secondary 
star may be the most reliable diagnostic we have of the 
white dwarf luminosity, since its intrinsic radiation is pro- 
duced in the far UV/soft X-ray regimes, where the effect 
of interstellar absorption will be very marked. If the irra- 
diating object is the white dwarf, the irradiation should 
decrease on the white dwarf cooling timescale. Prialnik 
(1986) shows how the surface layers of a white dwarf are 
heated during a nova explosion and cool as a power law on 
a time scale of 200 yrs. (This is in contrast to the cooling of 
white dwarfs after their initial formation, which involves 
cooling of the entire star and occurs on a time scale of 10 8 
years.) 

2. History of the photometric modulation in 
V1500 Cyg 

VI 500 Cyg was a naked eye nova, reaching a peak mag- 
nitude of V=2.2 in late August 1975. It is thought to be 
currently a slightly asynchronous AM Her system. The 
magnetic nature of the white dwarf was first proved by 
the detection of circularly polarized light (Stockman et 
al. 1988). The present day photometric period (orbital) 
is 1.8% longer than the polarimetric period (white dwarf 
spin). Presumably the two were thrown out of equality 
by the nova explosion and observations suggest that they 
will re-synchronize in a time scale of 200 years (Schmidt 
& Stockman 1991). 

Stockman et al. (1988) explain the photometric period 
evolution as follows: 

• The expansion due to the nova explosion increases 
the moment of inertia of the white dwarf, causing its spin 
period to increase to 0.141d, breaking the synchronous 
rotation. At this point the photometric modulation is as- 
sociated with the spin period of the white dwarf. 

• Interaction between the secondary star and the en- 
velope causes the envelope to be spun up to the binary 
period, strong coupling ensures that the core also achieves 
synchronism. 



M.W. Somers & T. Naylor: Post nova white dwarf cooling in V1500 Cygni 



• The remnant envelope shrinks back onto the white 
dwarf surface, reducing its moment of inertia and thus 
decreasing the spin period to 0.137d. 



Currently, therefore, V1500 Cyg displays two periods. 
There is a polarimetric signal at 0.137d which is the spin 
period of the white dwarf, and a photometric signal at 
0.140d which is the orbital period of the binary. Aside from 
Schmidt et al's spectrophotometry outlined in Section El 
further evidence that the secondary star now dominates 
the photometric modulation comes from the fact that the 
timing of flux maximum in our own data matches the or- 
bital ephemeris (see Section 0). It is also clear however, 
from the presence of nickering and the slightly asymmet- 
rical B band light-curves, that there is some contamina- 
tion of the orbital modulation. This could either be due to 
the accretion stream, or perhaps because the spin period 
of the white dwarf may also have a photometric signal 
(Pavlenko & Pelt 1988). 

3. Modelling white dwarf cooling after surface 
heating. 

Prialnik (1986) modelled the evolution of a classical nova 
through a complete cycle; accretion, outburst, mass loss, 
decline and resumed accretion. The model is for a 1.25M© 
C-0 white dwarf. The resulting outburst is fast and similar 
to that observed in V1500 Cyg, matching the composition 
of the ejected envelope very well. The white dwarf is mod- 
elled allowing for heat transfer via radiation, conduction 
(Iben 1975, Kovetz & Shaviv 1973) and convection (Mi- 
halas 1978). The result may be fitted well with a power 
law cooling curve of the form 



L oc t 



(1) 



where L is the luminosity of the white dwarf and t is the 
time since outburst. 

4. Observations 

To extend the baseline of photometric amplitude decline 
we obtained one orbital cycle of B and V band photometry 
using the JKT on La Palma on the night of 1995 October 
3. The TEK4 CCD was used with pixels binned 2 by 2 to 
achieve 0.6" xO. 6" pixels in rapid readout mode. The seeing 
was around 2.0 arcsec. The exposures were typically of 120 
seconds with filters being alternated between observations. 
A bias frame was subtracted from each image, and 
it was then flatfielded using an image of the twilight sky. 
The counts for V1500 Cyg and various other stars were ex- 
tracted from each frame using the optimal weighting pro- 
cedure described in Naylor (1998). We divided the counts 
for VI 500 Cyg by those for star CI of Kaluzny & Seme- 
niuk (1987), allowing us to put the lightcurves in Figures 
|l| and g onto a magnitude scale. 
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Fig. 1. The B band magnitude of V1500 Cyg plotted 
against Barycentric Dynamical Julian Date. 
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Fig. 2. The V band magnitude of V1500 Cyg plotted 
against Barycentric Dynamical Julian Date. 

The times of maximum and minimum agree, within 
errors, with the ephemeris of Semeniuk, Olech & Nalezyty 
(1995). The data do not show a pure heating modulation 
since there is some evidence of a dip in the light curve at 
time 2 449 994.61. Such dips are not uncommon, see for 
example the light-curves of Kaluzny & Semeniuk (1987) 
whose data are from 10 years earlier than our observations. 
These irregularities add errors in the estimation of the 
amplitude of the modulation (see Section ||) . 

5. The Amplitude v Time relationship 

By searching through various sources a record of the B 
band photometric behaviour observed in VI 500 Cyg since 
outburst has been assembled. Table 3.1 gives a list of all 
references used. The amplitude of the B band photomet- 
ric modulation versus the time since outburst is plotted 
in Figure on logarithmic scales. We immediately note 
that these points lie roughly in a straight line, implying 
a power law decay in time, although it should be noted 
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Table 1. The flux in the orbital modulation of VI 500 
Cyg. 



Days from 


Flux amplitude 


outburst 


(10~ 16 ergs cm~ 2 s _1 


320 


1060 ± 200 


329 


720 ± 190 


357 


310 ±32 


714 


35 ±6 


717 


26 ±6 


718 


30 ±8 


742 


21 ±5 


779 


39 ±6 


781 


38 ±8 


803 


34 ±6 


825 


24 ±8 


1072 


14 ±0.6 


1105 


12 ±0.6 


1172 


22 ±0.6 


1172 


18 ±0.6 


2167 


61 ±0.6 


2230 


12 ±0.6 


4039 


3.6 ±0.6 


4041 


2.9 ±0.3 


4305 


3.1 ±0.3 


4306 


3.9 ±0.3 


4307 


2.8 ±0.3 


4308 


3.1 ±0.3 


7339 


2.0 ±0.3 



A- 



Reference 

Kemp, Sykes & 

Rudy (1977) 
Patterson (1978) 



Home & 
Schneider (1989) f 

Kaluzny & 
Semeniuk (1987) 

Kaluzny & 
Chlebowski (1988) 



This Work 



^These points are from spectroscopy, band pass 4750-4790A. 



that the first year of data have been omitted. While other 
authors have attempted to find a trend in the amplitude 
versus time data they relied on the amplitude in magni- 
tude space. This is a measure of the amplitude relative to 
the brightness of other parts of the system. Clearly this is 
inappropriate for V1500 Cyg since the overall brightness 
in the B band changes due to a large variety of effects. The 
amplitude in flux space is a good probe of the irradiation 
since it does not need to be corrected for changes in rela- 
tive brightness of other parts of the system. Hence, what 
we refer to as the "flux amplitude" is obtained by con- 
verting the magnitudes at orbital maximum and minimum 
into fluxes (assuming a star of -B=0.0 gives 7.2 x 10 -9 ergs 
cm~ 2 s _1 A -1 ), and differencing them. An unweighted 
straight line to these data give 



Aozt 



-1.26±0.21 



(2) 



where A is the flux amplitude, and the error bar corre- 
sponds to la. 

It is apparent that there is an extra source of noise that 
acts on a short time scale compared with the general trend 
of amplitude decline. (Hence our choice of an unweighted 
fit.) This is likely to be due to the aforementioned orbital 
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Fig. 3. 
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The decaying photometric variability of VI 500 



dips shifting around the orbital light-curve. The position 
of maximum light from the accretion columns with respect 
to the secondary star's photometric hump will vary with 
beat phase. This will cause the apparently random excess 
scatter in the photometric amplitudes. 

6. The Irradiating Flux v Time relationship 

The bolometric flux amplitude from the secondary star 
will be proportional to the irradiating flux. However, the 
flux measured in any particular bandpass will not follow 
such a proportionality, as the bolometric correction will 
change with the temperature of the heated face. To extract 
how the irradiation is changing from the changing flux 
amplitude one must determine the response of the heated 
face of the secondary star to decreasing the heating. One 
may conveniently approximate the response as 



A oc Ff„ oc t XT >. 



(3) 



Here r\ is the power law decline predicted by Prialnik 
to be -1.14 while x is the response of the secondary star to 
heating. In fact x is not a constant but is a weak function 
of the temperature of the heated face of the secondary 
star. We therefore used the code described in Somers et 
al (1996) determine how the amplitude of the modulation 
changes with the irradiating flux. This model irradiates a 
Roche-lobe filling star from a point source at the position 
of the white dwarf. The irradiation then raises the tem- 
perature of each surface element in impinges on, such that 
all the incident energy is re-radiated. We will also assume 
that the radiation from the secondary star can be approx- 
imated as a blackbody, an assumption we shall examine 
later. We used a mass ratio q=3, inclination i=60° and 
underlying secondary star temperature of 3000K. These 
values are typical of those quoted for VI 500 Cyg and the 
results are, in any case, insensitive to the exact values. 

We began by irradiating a 3000K secondary star such 
that the flux when the irradiated face was towards the 
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observer was equivalent to a star of 8000K, matching 
the front face temperature observed by HST. This model 
should correspond to the last data point in Figure ||. The 
flux amplitude in Figure declines by a factor of 10 , and 
so we increased the irradiation in our model until the flux 
amplitude had increased by this factor. Over this range of 
interest, we found that 

A ex i£f (4) 

represented the data to better than 25 percent for all val- 
ues of Fi rr . We also used tried using the bolometric cor- 
rection and colours of model atmospheres given in Bessell 
et al (1998) instead of blackbodies to represent the flux. 
We found this changed x by less than 0.06 in the low irra- 
diation case, which is that most affected by the difference 
between model atmospheres and black bodies. 

With a value for x we can now use the observations to 
derive a value of -q by equating (g) to (g). This yields: 

?7 = 0.94 ±0.09. (5) 

This result is just (2.2c) consistent with the value Prial- 
nik (1986). find from purely theoretical considerations of 
r\ = 1.14. Especially given the nature of the approxima- 
tions made, this seems to support the conclusion that the 
photometric variation in V1500 Cyg, and by implication 
other old novae, is caused by irradiation from the white 
dwarf. 
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7. Conclusions 

The above shows that for at least the first 20 years from 
outburst the white dwarf cooling models match the avail- 
able observations. Given the Prialnik-type cooling law, 
with the observed value of 77, and the temperature of the 
irradiated face at some known time after outburst (from 
Schmidt et al. 1995) then we can calculate the typical time 
taken for irradiation of the surface to become negligible. 
The irradiation will drop off so that the incoming radia- 
tion is less than double the unheated surface luminosity of 
the secondary star about 280±140 years after the outburst 
of V1500 Cyg. Interestingly, we find that in WY Sge, now 
over 200 years since nova outburst, the irradiation from 
the white dwarf has declined to these levels (Somers et al 
1996), although in that system the disc is a complicating 
factor. Thus both V1500 Cyg, and WY Sge suggest that 
white dwarfs really do cool as the theory predicts. 
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